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Abstract 

Given a function / : X — > R on a topological space, we consider the preimages of intervals 
and their homology groups and show how to read the ranks of these groups from the ex- 
tended persistence diagram of /. In addition, we quantify the robustness of the homology 
classes under perturbations of / using well groups, and we show how to read the ranks of 
these groups from the same extended persistence diagram. The special case X = has 
ramifications in the fields of medical imaging and scientific visualization. 
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1 Introduction 

The work reported in this paper has two motivations, one theoretical and the other practical. 
The former is the recent introduction of well groups in the study of mappings / : X — > Y 
between topological spaces. Assuming a metric space of perturbations, we have such a 
group for each subspace A C Y, each bound r > on the magnitude of the perturbation, 
and each homological dimension p. These groups, and the diagrams that they generate, 
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extend the boolean concept of transversal! ty to a real-valued measure we refer to as robust- 
ness. Using this measure, we can quantify the robustness of a fixed point of a mapping ifTOl 
and prove the stability of the apparent contour of a mapping from an orientable 2-manifold 
to 191 ■ In this paper, we contribute to the general understanding of well groups by study- 
ing the real-valued case. Along the way, we also extend the general theory of well groups 
to incorporate relative well groups. Specifically, 

I. we give a general definition of relative well groups given a mapping / : X — > Y, a 
number r > 0, and a nested pair A' C A of subspaces of Y, and 

II. we characterize the relative well groups of / : X — !■ M whenever A is an interval and 
A' is a subset of the endpoints. 

Applications of this theoretical work are anticipated in medical imaging and scientific vi- 
sualization, where data in the form of real-valued functions is common. To mention one 
example, it is common to acquire information about internal organs through a magnetic 
resonance image, which results in a 3-dimensional array of intensity values, best viewed as 
a function from the unit cube to the real line. The predominant method for highlighting or 
extracting relevant substructures of this image uses preimages of real values. Generically, 
these are 2-manifolds, commonly referred to as contours or isosurfaces lfT2l . Sometimes, 
these 2-manifolds are complemented by preimages of intervals, referred to as interval vol- 
umes in visualization j l 11 . In this paper, we call the preimage of a value a level set, and 
the preimage of an interval an interlevel set, in which the interval can be closed, open, or 
half-open. We contribute to the state-of-the-art by 

III. explaining how the homology of level and interlevel sets can be read off the extended 
persistence diagram of the function, and 

IV. describing how the robustness of features in level and interlevel sets, quantified 
through well groups, can be read off the same diagram. 

Our results add up to a 'point calculus' in algebraic topology for mining the rich homolog- 
ical information contained in the extended persistence diagram of a real-valued function. 
The compactness of the data representation and the efficiency of the mining operations 
make the diagram an attractive graphical interface tool for studying 3-dimensional images. 
We view this tool as complementary to the contour spectra described in [I], which plot con- 
tinuously varying quantities, such as area and volume, across the family of level sets. The 
most novel aspect of our diagram is the robustness information, which has previously not 
been available. This novelty is combined with the unprecedented ease with which homo- 
logical information is accessible. There is also evidence for the practicality of the interface 
provided by the fast oct-tree implementation of the described concepts jj], which has been 
used to study 3-dimensional images of root systems of agricultural plants. 

Outline. In Section|2l we review necessary background on persistence, zigzag modules, 
and well groups. In Section[3] we explain the point calculus for interlevel sets. In Section 
m we extend the point calculus to include the robustness information provided by the well 
groups. Finally, Section|5]concludes the paper with a brief discussion of the contributions 
and of future research directions. 
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2 Background 



We divide the background material into three parts, introducing persistence and extended 
persistence in Section 12.11 explaining the extension to zigzag modules and level set pyra- 
mids in Section [Z!2l and defining absolute and relative well groups in Section l23] 

2.1 Forward Maps 

Traditional persistent homology is based on a nested sequence of spaces, which induces a 
linear sequence of homology groups connected by maps from left to right. We describe this 
concept in two steps. 

Persistence. The persistence of homology classes along a filtration of a topological space 
can be defined in a quite general context |f8l|. For this paper, we need only a particular 
type of filtration, one defined by the sublevel sets of a tame function. Given a real-valued 
function / on a compact topological space X, we consider the filtration of X via the sub- 
level sets Xr(/) = f~^{—oo,r], for all real values r. Whenever r < s, the inclusion 
Xr{f) ^ Xs(/) induces maps on the homology groups Hp(Xr(/)) Hp(Xs(/)), for 
each dimension p. Here we will use field coefficients so that the homology groups are 
torsion-free and are therefore vector spaces over the field. Often we will suppress the ho- 
mological dimension from our notation, writing H(X,.(/)) = 0^ Hp(X,.(/)); in this case, 
we will always assume that all maps H(Xr(/)) — > H(Xs(/)) decompose into the direct 
sum of maps on each factor A real value r is called a homological regular value of / if 
there exists e > such that the inclusion Xr-s{.f) ^ '^r+sif) induces an isomorphism 
between homology groups for every 5 < e. If r is not a homological regular value, then it 
is a homological critical value. We say that / is tame if it has finitely many homological 
critical values and if the homology groups of each sublevel set have finite rank. Assum- 
ing that / is tame, we enumerate its homological critical values j'l < r2 < . . . < r„. 
Choosing n + 1 homological regular values Si such that sq < ''i < < • ■ • < ''n < Sn, 
we put Xi = Xg. (/). The inclusions X^ ^ Xj induce maps f'-' : H(Xi) — )• H(Xj) for 
?^ and give the following filtration: 

= H(Xo) ^ H(Xi) ^ . . . ^ H(X„) = H(X). (1) 

We say a class a G H(Xi) is born at X^ if a ^ imf'~^'\ A class a born at X.; is said to 
die entering Xj if f'^'^a) £ imf'"^'-' but f''-'~-^(a) ^ imf*~^'-'~^. We remark that if a 
class a is born at X.;, then every class in the coset [a] = a + imf'^^'* is born at the same 
time. Of course, whenever such an a dies entering X^ , the entire coset [a] also dies with 
it. We represent [a] graphically as the point (r^, rj) in the plane. Drawing all birth-death 
pairs as points, we get diagrams like the ones sketched in Figures [T] and |3] Supposing that 
G M is different from all homological critical values, we collect all points in the upper- 
left quadrant defined by {b, b) to get all classes born before b and still alive; see the left 
diagram in Figure [1] Their number is the rank of the homology group of the sublevel set, 
rankH(Xfc(/)). 

Observe that we really need the extended plane to draw the points because some classes 
get born but never die, so the corresponding points have oo as their second coordinates. 
There is an elegant way around this minor annoyance, which we now describe. 
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Figure 1: From left to right: the ordinary, extended, and relative subdiagrams of Dgm(/). The 
number of points (not shown) in the dark shaded regions is equal to the rank of the homology group 
of the sublevel set defined by b. 



Extended persistence. Since the filtration in ([TJ begins with the zero group but ends 
with a potentially nonzero group, it is possible to have classes that are born but never die. 
We call these essential classes, as they represent the actual homology of the space X. To 
measure the persistence of the essential classes, we follow Q and extend the sequence in 
([T]i using relative homology groups. More precisely, we consider for each i the superlevel 
set X* = /^^[s„_i, oo). Note that we have X'^ = and X" = X by compactness. For 
i < j, the inclusion X* ^ X-* induces a map on relative homology H(X, X') — > H(X, X-'). 
These maps therefore give rise to the following extended filtration: 

0= H(Xo) ^ H(Xi) ^ ... ^ H(X„) = H(X,X°) ^ ... ^ H(X,X") =0. (2) 

We extend the notions of birth and death in the obvious way. Since this filtration begins 
and ends with the zero group, all classes eventually die. We also extend the graphical 
representation of the information contained by forming persistence diagrams, which we 
now introduce more formally. We have such a diagram for each dimension p; see Figure 
[T] Each diagram is a multiset of points in the plane, containing one point (r^, rj) for each 
coset of classes that is born at X^ or (X, X"^*+^), and dies entering X^ or (X, X"^-'+^). In 
some circumstances, it is convenient to add the points on the diagonal to the diagram, but in 
this paper, we will refrain from doing so. The persistence diagram contains three important 
subdiagrams, corresponding to three different combinations of birth and death location. 
The ordinary subdiagram, Ordp(/), represents classes that are born and die during the 
first half of (|2]i. The relative subdiagram, Relp(/), represents classes that are bom and die 
during the second half. Finally, the extended subdiagram, Extp(/), represents classes that 
are born during the first half and die during the second half of the extended filtration. Note 
that points in Ordp(/) all lie above the main diagonal while points in Relp(/) all lie below. 
On the other hand, Extp(/) may contain points on either side of the main diagonal. By 
Dgm(/), we mean the points of all diagrams in all dimensions. Drawing these subdiagrams 
side by side can be cumbersome, and drawing them on top of each other can be confusing. 
In Section[3] we will introduce a new design that addresses these concerns. 



4 



2.2 Mixed Maps 



We note that the homology groups in the extended filtration of (|2|i, or in the shorter filtration 
of ([T]i, are all vector spaces over a fixed field and that the maps between them are all linear 
maps. In jSj, Carlsson and de Silva generalize this situation to sequences of vector spaces 
that are connected by maps going from left to right or from right to left. We now briefly 
review their work as well as the related work on level set zigzag modules in 16]. 

Zigzag modules. A zigzag module W is a finite sequence of vector spaces connected by 
linear maps which either go forward or backward between consecutive spaces: 

Wi ^ W2 o . . . o o Wj+i ^ . . . o W„. (3) 

If the arrow advances from Wj to Wj+i, then we denote the corresponding linear map as 
aj : \Nj — > Wj+i; otherwise, we write bj : W^+i — > \Nj. A submodule U of W is a 
collection of linear subspaces Uj C \Nj such that aj(Uj) C Uj+i or &j(Uj+i) C Uj, 
whichever is the case for j. A submodule U is a summand if there is a complementary 
submodule V, meaning every vector space splits as a direct sum = Uj ® V^. The 
authors in ||5| prove that every zigzag module can be split into indecomposable summands 
of a certain form, and, in particular, it has a basis, a concept we now describe. First, 
we suppose that we have, for each j, a set of elements £ Wj such that the nonzero 
elements form a basis of Wj . In other words, we can decompose Wj into the direct sum 
Wj — noting that some of the terms on the right hand side may be zero. We 

use the superscripts to form correspondences between the bases. Specifically, we require 
aj(Uj) = u^jj^i, or = u^, depending on the case. Furthermore, we assume that, 

for each superscript i, there exist x < y such that ^ iff j G [x,y]. In other words, for 
each fixed i, we have a submodule 

{u\) ^ (u^) ^ ...^ {u)) ^ ^ ...^ (u^) (4) 

of W in which the non-zero vector spaces are 1 -dimensional and form a single contiguous 
subsequence connected by identity maps. Calling such a submodule an interval module, we 
think of it as being in correspondence with the closed interval [a;, y\. The collection {u}} 
is a basis for the zigzag module if W can be decomposed into the direct sum of the interval 
modules Equivalently, the collection is a basis for W if each map aj is the direct sum 
of the maps (u*) {u^j+i), and each map hj is the direct sum of the maps — > (w*), 
whichever one is defined. 

Although a zigzag module W can have many different bases, the set of intervals asso- 
ciated to any such basis will be unique Q. For example, any basis for the zigzag module 
given by the filtration in ^ will have one interval [x, y] for each coset of classes bom at 
Xa; and dying entering \y. 

Mayer- Vietoris diamonds. We are interested in an elementary operation that connects 
two minimally different zigzag modules: a Mayer-Vietoris diamond. We suppose that we 
have two zigzag modules differing only at position j, and that at this position we have a 
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diamond of the following form; 



H(¥,V') 




H(C,C') H(D,D') ^ , (5) 




H(E,E') 

where we show the more general, relative form in which the primed spaces are subspaces 
of the corresponding unprimed ones, and we have E = C H B, E' = C fl D', V = C U B, 
and V = C U B'. We get the more special, absolute form by setting C = B' = E' = 
V = 0. The name of the diamond is justified by the long exact sequence we get by reading 
the diamond from bottom to top and iterating through the dimensions. When the primed 
spaces are all empty, this gives the classic version of the Mayer- Vietoris sequence, and 
more generally, we get the relative version: 

... ^ Hp(E,E') ^ Hp(C,C')©Hp(B,D') ^ Hj,(V,V') ^ Hp_i(E,E') ^ 

see e.g. |fT3l . Importantly, this sequence is exact, which means that the image of each map 
equals the kernel of the next map. 

Such diamonds arise in the following context. Consider again the function / : X — > M 
and the interleaved sequence of homological regular and critical values: sq < '"i < si < 
... < r„ < s„. Setting H{f-\sj)) and Wsj+i = H(/-i[sj, s^+i]), we get a 

zigzag module of length 2n + 1, which, following 0, we refer to as the level set zigzag of 
/. It starts and ends with and alternates between advancing maps a2j and backward 
maps &2j+i- From this module, we can create a new one by fixing an index j, sub- 
stituting [sj,Sj+2] = [sj,Sj+i] U [sj+i,Sj+2] for Sj+i = [sj,Sj+i] n [sj+i,.Sj+2], and 
leaving all other groups unchanged; of course we also reverse the two maps involving the 
changed space. This produces a new zigzag module which differs from the old via a Mayer- 
Vietoris diamond. This construction can be generalized by flipping between intersections 
and unions of larger intervals and pairs of intervals, thus producing a whole array of zigzag 
modules which differ via Mayer- Vietoris diamonds. 

The pyramid. Starting with the level set zigzag, we get an array of zigzag modules which 
are best described as monotonic paths that go diagonally up and down, always from left to 
right. The array of such paths is connected within a pyramidal structure, which we now 
describe. As a graphical guide, we consider the square drawn in Figure |2] We give it a 
coordinate system by parameterizing the downward slope from co at the upper left, to — cx) 
in the middle, and back up to oo at the lower right. Similarly, we parameterize the upward 
slope from — cxo at the lower left, to oo in the middle, and back to — oo at the upper right. 
The two slopes divide the square into four triangular regions, each containing a point with 
coordinates a and b for every choice of a < 6. We interpret this point differently in each 
of the regions. To explain this interpretation, it is convenient to introduce a shorthand that 
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uses open set notation for pairs of closed sets, writing A — A' for (A, A'). Specifically, 

r\x,y) = (/-i(-^,^),/-i(-oo,x] U/-i[2/,^)). 
If a point with coordinates x and y lies in the bottom region, we think of it as the space 

z y X w 




w X y z 

Figure 2: Points in the pyramid are absolute and relative homology groups. Monotonia paths are 
zigzag modules, any two of which differ by a finite number of Mayer- Vietoris diamonds. 

f~^[x,y]. However, if the point lies in the left, right, or top region, we think of it as 
f~^{x,y], f~^[x,y), or f~^{x,y), respectively. If we now take w < x < y < z and 
consider the points (lu, y), (w, z), (x, y), and (x, z), we get a Mayer- Vietoris diamond in 
each region; see Figure |2] This is easiest to see in the closed interval case since [x^y] = 
\w,y\ n \x,z\ and \w,z\ = [w,y\VJ[x,z\. In the closed-open case, we have [x, oo) = 
[w, oo) n [x, oo) and [w, 00) = [w, 00) U [x, 00) as well as [z, 00) = [z, 00) n [y, 00) and 
[y, 00) = [z, 00) U [y, 00). Similar computations verify the diamond in the remaining two 
cases. 

By repeated application of the diamond, we can generate any monotonic path from the 
one along the bottom edge of the square. Each path is thus decorated by spaces as de- 
scribed, and applying the homology functor gives a zigzag module of absolute and relative 
homology groups. The latter arise when we move the left or right end of the path, which 
can be done without the Mayer- Vietoris diamond because the corresponding spaces are and 
stay empty so that the module remains unchanged. Besides the level set zigzag along the 
bottom edge, we are particularly interested in the path along the upward slope, which trans- 
lates into the extended filtration of Its midpoint is (—00, 00), the center of the square. 
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which results in H(/ ^{—oo, oo)) = H(X). For this reason, we think of the center as the 
apex of a pyramid, as viewed from above. 

Remark. As a partial justification for the notation with open sets, we mention that the 
homology group of the preimage of the interval {x, y), if computed with infinite chains, is 
isomorphic to the relative homology group of (/^^ [x, y] , f^^{x) U f^^{y)). By excision, 
this is isomorphic to the relative homology group of (/^^(— oo, oo], /^^(— oo, a;] U f~^[y, oo 

2.3 Perturbations 

The reader who wishes to learn how to read the homology of interlevel sets can safely skip 
Section 12.31 and now continue with Section |3] However, to differentiate the robust from 
the non-robust homological information in these readings, we need to first understand the 
subgroups of homology that give meaning to this concept. 

Well groups. Suppose that we have a continuous mapping / : X — > Y between topo- 
logical spaces. Given a subset A C Y, we review here the definition of the well groups 
Ua(?') for each radius r > 0. When A is clear from context, we will drop it from the 
notation and simply write U(r), by which we mean the direct sum of groups Up(r), for 
each homological dimension p. We will also need the assumption that /~^(A) has homol- 
ogy groups of finite rank in each dimension. In addition to the mapping /, we assume a 
subspace V of C(X, Y), the space of continuous mappings from X to Y, requiring that V 
contains /. For example, V might consist of all mappings homotopic to /. We assume 
a metric on V and write ||/ — h\\^ for the distance between two mappings. We call h an 
r-perturbation of / if ||/ — < r. Given A C Y, we introduce the radius function, 
/a : X — R, by setting fhix) to the infimum value of r for which there exists an r- 
perturbation h G V with h{x) G A. We filter X via the sublevel sets of the radius function, 
setting X, (/a) = ./a'[0,'^]. For r < s, there is a map f^-^ : H(X,(/a)) ^ H(X4/a)). 
The preimage of A under any ?- -perturbation h of f will obviously be a subset of Xr(/A), 
and hence there is a map on homology, j/i : H(/i~^(A)) — > H(Xr(/A))- Given a class 
a € H(Xr(/A)) and an r-perturbation h of /, we say that a is supportedhy /i if a G im.jft,. 
The well group U(r) C H(Xr(/A)) is then defined ifTOl to consist of the classes that are 
supported by all r-perturbations of /: 

U(r) = fl imj;,. 

Il''-/llp<'- 

For r < s, the map f^''* restricts to a map U(r) — > H(Xs(/a)). On the other hand, 
H(Xs(/a)) contains U(s) as a subgroup. It can be shown that U(s) C f]^'^(U(r)) whenever 
r < s; see |[TOl . In other words, the rank of the well group can only decrease as the thresh- 
old value increases. We call a value of r at which the rank of the well group decreases 
a terminal critical value of /a. The well diagram of / and A is the multiset of terminal 
critical values of /a, taking a value k times if the rank of the well group drops by k at the 
value. Often we will refer to this diagram as the robustness of the preimage /"^(A). In 
this paper, we focus on the case Y = M and V — C(X, M), lifting the usual metric on 
E to P by defining ||/ — h\\^ = \\f — h\\^ = sup^g^ \f{x) ~ h{x)\. In this case, the 
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radius function satisfies /A(a;) = infaeA 1/(2;) ~ a|- In general, the relationship between 
the terminal critical values and the homological critical values of /a is not completely un- 
derstood. However, if Y = M and A is a point, we will see that the former is a subset of the 
latter We get more complicated relationships when A is an interval. 

Example. Consider the torus X, as shown in Figure[3] along with the vertical height func- 
tion / : X — > M and the space A = {a}. The preimage of A, /^^(A) = ./^ ^(0), consists 
of two disjoint circles on the torus; hence there are two components and two independent 
1-cycles, all belonging to the well group at radius 0. For small values of r, Xr(/A) consists 
of two disjoint cylinders. The homology has yet to change; furthermore, although the proof 
will come later, all classes still belong to the well groups at these small radii. 




Figure3: Left: the torus and the preimage of the interval [a~r,a+r]. Right: the extended persistence 
diagram of the vertical height function. Each point is labeled by the dimension of the corresponding 
homology class. The dark shaded portions of the diagram represent the homology of /^^ [a — a+r]. 

Now consider the value of r shown in Figure |3] For this r, the sublevel set X^ = 
Xr(/A) consists of two pair-of-pants glued together along two common circles. We note 
that Ho(Xr) has dropped in rank by one, while the rank of Hi(Xr) has grown to three. In 
contrast, the rank of Ui(r) is less than or equal to one. Indeed, the function h : X ^ W, 
defined by /i = / — r, is an r-perturbation of / and the zero set of the corresponding 
distance function, /i^^(O) = f~^{a + r), is a single closed curve. Since the rank of the 
first homology group of that curve is one, and since the rank of m\\h can be no bigger than 
this rank, the well group U i (r) can also have rank at most one. That it does in fact have 
rank exactly one will follow from our results in Section]?] 

Relative well groups. Since the pyramid involves relative homology groups, it seems 
wise to extend the definition of well groups into the context of relative homology. While 
this notion is new, it follows the above ideas closely so that presenting the definition in this 
background section seems appropriate. Assume again that we have a continuous mapping 
/ : X — > Y between topological spaces, as well as a subspace V of C(X, Y) that contains / 
and is equipped with a metric. Given a nested pair A' C A of subspaces of Y, and a radius 
r > 0, we note that X'. = X,(/a') is ^ subset of X^ =Xr{fA)- For each j'-perturbation /i 
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of /, there is an inclusion of pairs (/i^^(A), h^^{A')) ^ (X^, XJ,), which induces a map 
jh '■ H(/i~^(A), /i~^(A')) — > H(Xr, XJ,) between relative homology groups. The relative 
well group U(A,A')('') defined to be the intersection of the images of these maps, taken 
over all r-perturbations of /: 

U(A,A')('') = n 

ll'»-/llp<'- 

When a distinction is needed, we will refer to the previous notion of well groups as absolute 
well groups. 

3 Combinatorics of Homology 

In this section, we present the first half of our point calculus, showing how to read the 
homology of a level or interlevel set from the extended persistence diagram. The crucial 
technical concept is that of a basis of the pyramid of zigzag modules, which we establish 
by strengthening the Pyramid Theorem in ||6] . 

Flipping a basis. We construct a basis for the pyramid one step at a time, by flipping the 
basis of one zigzag module to the next. For this purpose, we consider two zigzag modules 
that differ at one position, and we assume that there is a Mayer- Vietoris diamond serving 
as a connecting bridge between the two modules at that position. Drawing the diamond 
with the intersection at the bottom and the union at the top, as in Q, we say the diamond 
connects the lower module with the upper module. Given a basis of the lower module, we 
can show that we can construct a basis of the upper module so that the two bases agree on 
the overlap. We refer to this operation flipping the first basis to the second. 

Basis Flip Lemma. Given two zigzag modules that differ by a single Mayer- Vietoris 
diamond, we can flip any basis of the lower module to a basis of the upper module. 

Proof. We give a proof by construction. Writing {e^} for the basis of the lower zigzag 
module, we describe a basis of the upper zigzag module that differs from the lower 
one only at the position j at which the modules differ; as in (|5]l. We thus at once set v], = e\ 
for all k 7^ j, and the main task is then the construction of the u*. Put briefly, our rule will 
be that ^ iff an odd number of e*_]^, ej, e*^j^ are non-zero. We give more specifics 
via a case analysis. The cases are labelled pictorially, with black dots denoting non-zero 
classes, showing only the positions j — 1, j, j + 1- 

Case 1 (• -^-^ )'■ We have ej_]^ ^ and — e'j^i ~ 0, and define as well as the 
advancing map using the Mayer- Vietoris diamond, namely = aj-i{e^j_i), which 
is non-zero by exactness and because e* = 0. 

Case 2 (k^ ): Again we set Vj ~ aj_i(e*_j), which is zero by exactness and be- 
cause e* 7^ 0. 

Case 3 ( v^^^): We set = aj-i{ej_i) = which in this case is non-zero. 

Indeed, if it were zero, then, by exactness, the pair {e)-i, 0) would be in the image 
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of the map 6j_i © aj. By the direct-sum decomposition of the maps in the basis, this 
would imply that aj{ej) = 0, a contradiction. 

Case 4 (_->'): We have ^ and e^„i = e*^^ = 0. If there are £ > indices i of 
this kind, then the orthogonal complement to the image of the map cj, defined below, 
has rank £, as we prove shortly. We pick £ classes that span this complement. Since 
Vj maps to e* via the connecting homomorphism of the Mayer- Vietoris sequence, the 
homological dimension of is one higher than that of e* . 

Case 5 ( ^ ^ ■): This is symmetric to Case 2, and we set — bj{e^j^i) ~ 0. 

Case 6 ( • ^ This is symmetric to Case 1, and we set — foj(ej+i) ^ 0. 

Note first that we now have interval modules {e^ } in the lower zigzag module, and interval 
modules {vj} in the upper zigzag module. To show that the latter are indeed summands, 
we only need to verify that the non-zero classes Vj form a basis of H (V, V), the new group 
in the upper zigzag module. Using the notation in (|5]l, we let E denote the vector space 
spanned by the pairs {e'j_l,e'j_^_l), noting that E is a subspace of H(C, C) © H(D, W), but 
because of Case 3 it is not necessarily the entire direct sum. We consider the subspaces E jv 
of E spanned by the pairs {e'j_i, e]+i) in each Case N, for 1 < iV < 6. These subspaces 
are independent and span the entire space E. In other words, zero is the only element 
common to any two of the subspaces, and the ranks of the subspaces add up to the rank of 
E. 

The case analysis suggests a map cj : E ^ H(V, V) with Cj((e*_]^, e]_|_]^)) = w], if 
(e*_]^, e*_|_j) ^ (0,0), and zero otherwise. Since E4 — 0, this map is zero on E4, but it 
is also zero on E2 and E5. Furthermore, cj is injective when restricted to Ei, E3, and Eg. 
We proceed to show that the images of these latter three vector spaces under Cj are inde- 
pendent of one another To derive a contradiction, we first suppose that Cj(Ei) n cjIEq) 
contains a non-zero class. Then there must exist {a, 0) S Ei and (0, /3) G Eg with 
aj-i{a) = bj{/3) ^ 0. Hence, {a, 13) G ker(aj_i (Bbj), which, by exactness, tells us 
that a G imbj-i. But this contradicts the direct-sum decomposition of the map 
Next, suppose that Cj(Ei) n Cj(E3) contains a non-zero class, which means there exists 
(a,0) G El and (7,/?) G E3 such that aj-i{a) = bj{f3) 7^ 0. As above, this implies 
that (a, /?) G ker (aj_i © bj), and we reach the same contradiction. Finally, a symmet- 
ric argument gives Cj(E3) n Cj(Eg) = 0. We conclude that Cj(Ei), Cj(E3), and Cj(Eg) 
are independent subspaces of H(V, V). In Case 4, we picked a basis for the orthogonal 
complement to their span; all together, we have a basis of H (V, V), as required. H 



Establishing a basis. The Pyramid Theorem in ||6] establishes an explicit bijection be- 
tween the interval modules that arise in the decomposition of any two zigzags within the 
pyramid. We strengthen this result by establishing bases on all the zigzag modules in such a 
way that the basis elements correspond to the intervals and respect the same bijections. We 
call this a basis of the pyramid. To construct it, we note that the paths in the pyramid are 
connected by Mayer- Vietoris diamonds. We can therefore flip a basis of the level set zigzag 
upwards through the entire pyramid via repeated application of the Basis Fhp Lemma. 
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Pyramid Basis Theorem. A basis of the level set zigzag module extends to a basis 
of the entire pyramid. 

We now give an explicit description of how the interval modules of the various paths in 
the pyramid relate to each other A convenient reference in this description is the extended 
filtration which follows the upward slope through the middle of the pyramid. Its first 
half is parameterized from — oo to oo, and its second half from oo back to -co. Let now 
X and y be two points along the upward slope, with x to the left of y. We distinguish 
between the ordinary case (x < y, both in the first half), the relative case (y < x, both in 
the second half), and the two extended cases [x < y and y < x, with x in the first half 
and y in the second half). For each case, we sketch how the basis element of the interval 
corresponds to basis elements of other homology groups in Figured As a general pattern, 




extended relative 

Figure 4: The basis element that corresponds to the interval from x to y along the upward slope maps 
to all spaces between the paths of its two endpoints. The four squares show the pattern for the four 
different types of intervals. 

the two points trace out two curves consisting of segments with slopes ±45° that reflect 
before they hit the vertical sides and end at the horizontal sides of the square. The reason 
for the slopes are Cases 1, 2, 5, and 6 in the proof of the Basis Flip Lemma, and the reason 
for the reflection is the local change in the zigzag structure caused by moving the terminal 
zero group up. The two curves cross at one point inside the square, and the location of 
that point is characteristic for the case (the triangular region on the left in the ordinary 
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case, at the top and at the bottom in the two extended cases, and on the right in the relative 
case). The crossing is caused by Case 4, in which the correspondence between the basis 
elements is constructed via the connecting homomorphism of the Mayer- Vietoris sequence 
and therefore comes with a shift in homological dimension. 

Turning the table. The regions in Figure |4] show all the spaces represented by points in 
the pyramid to which the basis element corresponding to the interval [a;, y] is relevant. We 
are now interested in the inverse question: which basis elements are relevant to a given 
space? More specifically: which intervals in the decomposition of the extended filtration 
(|2]i map to the basis of the homology group of the space represented by a point with coor- 
dinates a and 6? We answer this question by considering the following subregions of the 
p-dimensional persistence diagram: 

Xp[a,b] = {(.T,y) e Ordp(/) I a; < 6 < y} U {(a;,y) e Extp(/) I a; < 6,a < y}, 

Qp[a,b] = G Extp(/) I 6 < a::,y < a} U {(a;, J/) e Rclp(/) I y < a < x}, 

Ap[a,&) = G Extp(/) I a < y < 6} U {(x-,?/) e Relp(/) I a < y < 6 < 4, 

Qp[a,b) = {(.T,y) e Relp(/) I y < a < a: < 6}, 

Ap(a, 6] = {(a;,y) e Ordp(/) | a; < a < y < 6}, 

gp(a,6] = {(.T,y) G Ordp(/) I a < a; < 6 < y}U {(a;,y) G Extp(/) I a < a; < 6}, 

Ap(a,&) = {(x,y) G Ordp(/) I a; < a < y} U {(a-,y) G Extp(/) I a; < a,fe < y}, 

£»p(a,&) = {(x,y) G Extp(/) I a < a;,y < fe} U {(a;,y) G Rclp(/) I y < 6 < x}, 

where we assume that a and b are both homological regular values. These multisets are 
displayed in Figure |5] in which we have also introduced a new, and for our purposes more 
convenient, way of drawing the extended persistence diagram. We have glued the domains 
of the three sub-diagrams and drawn the result as a right-angled triangle. In this triangle, 
the birth and death axes go from — oo up to +oo and then continue on back to — oo. In other 
words, we flip the extended subdiagram upside down and glue its (formerly) upper side to 
the upper side of the ordinary subdiagram. Similarly, we rotate the relative subdiagram 
by 180 degrees and glue its (formerly) right side to the right side of the (flipped) extended 
subdiagram. After gluing the three domains, we rotate the design by —45 degrees so the 
triangle rests on its longest side, consisting of the diagonals in the ordinary and relative 
subdiagrams. The diagonal of the extended subdiagram is now the vertical symmetry axis 
passing through the middle of the triangle. 

Remark. There is a straightforward translation of this triangular design to the repre- 
sentation of persistence advocated in 14]. Namely, draw a isosceles right-angled triangle 
downward from each point in the multiset and call the horizontal lower edge the corre- 
sponding bar. The barcode is the multiset of bars, one for each point in the diagram. 
Similarly, we can translate the triangular design into the square design of the pyramid by 
cutting along the vertical axis, turning the right triangle upside-down, and gluing the two 
triangles along their hypotenuses. 

Reading interlevel sets. The purpose of the multisets defined above is to offer a conve- 
nient way to read the absolute or relative homology of an interlevel set from the extended 
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H(rHa,61) 



H(.r'[a,6)) 




H(.rH«,61) 



H(.r'(a,6)) 



Figure 5: The triangle design of the persistence diagram showing the regions A and g for the four 
types of intervals in darker shading. When we collect the points to compute the rank of the p-th 
homology group, we shift the homological dimension of classes as shown. 

persistence diagram. We need some definitions to combine all four types into one. First, 
we let B be the collection of interval modules in the decomposition of the extended filtra- 
tion (O. As mentioned earlier, this collection is in bijective correspondence with the points 
in Dgm(/). We write V = {B) for the abstract vector space spanned by B, and we let 
V = {{B') \ B' Q B} be the collection of vector spaces spanned by subsets of this basis. 
Second, we write 



for the region of points in the persistence diagram that correspond to the basis elements of 
Hp(/~^ (/)), and call it a pair of wings. With these concepts, we have the following result, 
which implies that the rank of Hp(/~^(/)) is the number of points in Wp(/): 

iNTERLEVEL SET LEMMA. For each dimension p and each interval / whose endpoints 
are homological regular values, there exists an isomorphism that takes Hp(/^^(/)) to the 
vector space Gp (/) G V spanned by the basis vectors corresponding to the points in Wp (/) . 

Proof. Write B = {e'} and let {w*} be the basis of the group Hp(/"^(/)), where / is an 
interval with endpoints a < b that can be closed, closed-open, open-closed, or open. The 
claimed isomorphism is then the linear map 7 : Hp(/~^(/)) V defined by 7(w*) = {e*} 
for all non-zero v\ 




Xp[a,b]U gp+i[a,b] if I 

Xp[a,b) U gp+i[a,b) if I 

Xp^i{a,b] U gp{a,b] if I 

Xp^i{a,b)U gp{a,b) if I 



[a,b], 
[a, 6), 
(a, 6], 
(a, 6), 
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To understand why the image of 7 consists of the intervals that correspond to the points 
in Wp(/), we need to recall the transformation rules sketched in Figure |4] Consider for 
example the closed interval case, / = [a, b], for which Wp{I) — Xp [a, b] U Qp+i [a, b]. Since 
the interval is closed, the homology group is represented by the point (a, b) in the lower 
triangular region. To lie in the dark shaded region, this point must satisfy the constraint 
X < b < y in the ordinary case, x < b and a < y in the first extended case, and x < b 
and a < y without dimension shift in the second extended case. These inequalities define 
Xp[a, b]. Furthermore, we get b < x and y < a with dimension shift in the second extended 
case, and y < a < x, again with dimension shift, in the relative case. These inequalities 
define Qp+i [a, b], which completes the proof in the closed case. For a proof of the closed- 
open, open-closed, and open cases, note that the points representing Hp(/~^(/)) are found 
in the right, left, and top triangular region of the pyramid, and then argue in a similar 
fashion. H 



4 Combinatorics of Robustness 

The definition of well group given in Section|2]involves an uncountable number of pertur- 
bations, which give rise to the intersection of a potentially large number of subgroups, and 
as such does not seem amenable to computation. In this section, we show that the situation 
in the real-valued case is simpler, and that we are able to read the absolute and relative well 
groups directly from the extended persistence diagram. We begin with a consequence of 
the Mayer- Vietoris sequence, which provides the main technical ingredient of our proofs. 

A corollary of Mayer- Vietoris. For convenience, we establish the following notational 
convention, wherein we reuse the same letter in different fonts. We will need it for abso- 
lute and for relative homology groups. To avoid repetition, we state it now for the more 
general relative case. Letting U' C U and V C V be pairs of topological spaces, we 
write (U,U') ^ (V,V') if U C V and U' C V'. This inclusion of pairs induces a map 
u : H(1LJ, U') — > H(V, V) on homology groups, and we write U = imu for the image of 
this map. Note that U is always a subgroup of H(V, V), namely the subgroup of homology 
classes that have a chain representative carried by (U, U'). Note also that the rank of U can 
never exceed the rank of H(V, V). Suppose that, furthermore, (T, T') ^ (U, U'). Then, 
from the sequence of maps H(T,T') H(U, U') ^ H(V, V), we see that T, the image 
of H(T, T') in H(V, V), must be a subgroup of U. The following lemma is a direct con- 
sequence of the exactness of the Mayer- Vietoris sequence. However, we will use it often 
enough that it seems reasonable to state and prove it formally. 

Mayer-Vietoris Lemma. Suppose the pair of topological spaces V C V can be 
decomposed as V = C U © and V = C U where C C C and D' C B. Set (E, E') = 
(C n D,C' n ©')■ If a class a e H(V, V) belongs to C and to D, then a also belongs to E. 

Proof. Following our convention, we use the notation c : H(C, C) — > H(V, ¥') for the 
map on homology induced by the inclusion of (C, C) in (V, V). Similarly, we write d : 

H(©,D') ^ H(V,V') ande : H(E,E') ^ H(V,V'), as well as : H(E,E') ^ H(C,C') 
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and Bd : H(E,E') ^ H(©,P'). Note that C = inic, D = imd, and E = ime. Consider 
now the relevant portion of the Mayer- Vietoris sequence for (V, V'): 

H(E,E') ^-^^ H(C,C') © H(D,P') H(V,V'). 

By assumption, a G C, so there exists some ac £ H(C, C) such that c(ac) = a. Similarly, 
there exists an ad G H(D,D') such that d{ad) = ol. This implies that the pair (ac,ad) 
belongs to the kernel of c — d, and thus also, by exactness of the sequence, belongs to the 
image of (Sc, e^). Hence, there exists a,; G H(E,E') with ^cip^e) ~ etc and ^dipLe) — cud- 
In particular, since e = c o Sc, we have e{ae) ~ a, and therefore a G E as claimed. H 

In the typical application of the Mayer- Vietoris Lemma, we will construct further pairs 
(T, T') ^ (C, C) and (B, 1') (D, ©') such that a G T n B. From the remark above, 
we know that T C C and BCD. The lemma then applies and we can conclude that a G E, 
as before. 

The well group of a level set. As a warm-up exercise, we first consider the case in which 
A is a single point. More specifically, we suppose that we have a compact topological space 
X and a function / : X — 5- R, and we find the well groups U(r) = Ua(?'), where A = {a} 
is some point on the real line. In this case, X,.(/a) = /a ^[0; ^] = f~^[o- — r,a + r]. To 
state the formula, we distinguish two particular subspaces of X^ = X,.(/a), namely the 
top level set, T,- = f^^{a + r), and the bottom level set, B,. = f^^{a — r). Using the 
convention from before, we write and B^ for the images of H(Tr) and H(Br) in H(Xr). 

Point Formula. U(r) =7^0 Br, for every r > 0. 

Proof. We prove equality by establishing the two inclusions in turn. To show U(r) C 
Tr n Br, consider an arbitrary class a G U(r). We define /itop ^ f — r and ft-bot = f + r 
and note that they are r-perturbations of /, with ft.^p(a) = Tr and ^■bo\(a) = Br. By 
definition of the well group, a is supported by every r-perturbation of /, and therefore by 
/itop and by /ibot- It follows that a G Tr fl Br. To show Tr fl Br C U(r), we consider an 
arbitrary class a G Tr n Br and let h be an arbitrary /--perturbation of /. To finish the proof, 
we need to show that a is supported by h. We define Cr — h^^[a, cx)) fl X.,. and Br = 
/i"i(-oo, a] n Xr. Note that Cr U Br = Xr while Cr n Br = /i"^(a). Furthermore, the 
inequality — /Ho^ < r implies that Tr C Cr and Br C Br. By the Mayer- Vietoris 
Lemma, a is supported by h^^{a), as required. H 

Remark. The Point Formula implies that the well group for a Morse function / can 
change only at critical values of the function /a, where A = {a}. In other words, terminal 
critical values are, in this simple context, just ordinary critical values. Indeed, if [r, s] is 
an interval that contains no critical values of /a, then there is a deformation retraction 
Xs(/a) Xr(/A) providing an isomorphism H(Xs(/a)) H(Xr(/A))- Furthermore, 
this retraction maps Tg onto Tr, in such a way that that the images of H(Tr) and H(Ts) 
in H(Xs(/a)) are identical. Similarly, the images of H(Br) and H(Bs) in H(Xs(/a)) are 
identical. Hence the well groups U(r) and U(s) are isomorphic. 
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The well group of an interlevel set. We generalize from a point to an interval, which 
can be closed, closed-open, open-closed, or open. To that end, we define the spaces and 
maps so that the formula for the well group is the same in all four cases, and indeed the 
same as in the Point Formula above. Assume a < 5, set A = [a, b], and let A' C {a, b}. 
We thus getXr = Xr(/A) = f~^[a — r,b + t] and ~ X^ ( /a' ) > which is the empty set, 
f~^[b—r, b+r], f~^[a—r, a+r], or the union of these two interlevel sets. Correspondingly, 
we define the top and bottom interlevel sets: 

T, = f-^[a + r,b + r], T; C {J-\a + t), J-^b + r)}, 

1, = r>-r,b-r], B; C {r\a-r)J-\b-r)}- 

see Figure|6] The pairs (T,., T^) and (B^, B^) include into (X^, X^) in all four cases. Still 



1 1 


1 1 } 1 










1 




1 




h+r 
h 

b-r 


















1 




1 




a + r 
a 

a — T 






1 — ' — r 1 1 


1 



Figure 6: Each vertical strip represents X, and the shaded portions mark (Cr, Cj.) and (Tr, T^) on 
the left, (X,., X^) in the middle, and (B,., B^) and (D^, D^) on the right. 

using the notational convention from above, we write T^ and for the images of H (T^ , ) 
and H(Br,B;) in H(X^,X;). The formula for the well group, U(r) = U(A,A')(f), is then, 
unsurprisingly: 

Interval Formula. U(r) = 7^0 B^, for every r > 0. 

Proof. We give the argument for the most complicated of the four cases, when A' = 
{a, 6}. The proofs of the other three cases are simpler versions of the same argument. 
We may assume a + r < b — r, else X^ = X^., which implies that all groups in the 
claimed formula are zero and so we are done. To prove the inclusion U(r) C T^ fl B^, we 
consider the two r-perturbations /itop ~ f — r and /ibot = f + r, as before. Note that 
(Tr, T^) ~ h^op{a, b) and (Br,B^) = ''■bot('^' desired inclusion follows from 

the definition of relative well groups. To prove T^ n B^ C U(r), we choose an arbitrary 
class a e Tr n Br and an r-perturbation h of /. Furthermore, we introduce the following 
pairs of subspaces: 

Cr = /i"^[a, oo) n /~^(-oo, 6 + r], 

C; = {h-'^[a,oo) n f-\-oo,a + r]) U {h-'^[b,oo) n f-\-oo,b + r]), 

Br = h-^{-oo,b] n f-^[a-r, oo), 

B; = ih-'^{-oo,a]r\ f-'^[a-r,oo)) U {h-\-oo,b]r\ f-^[b-r,oo)); 
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see Figure |6] for a depiction of the open case. Since h is an r-perturbation, we have 
(Tr,T;) ^ (Cr,C;) and similarly (B^,B;) ^ (Dr,©^)- This implies T,. C and 
Br C Dr, and therefore a e n D^. It is easy to see that (C^ U Dr,CJ. U D^) = 
(X^,X;), and also that (C^ n Br,C'r n = {h~^{A),h~^{A')). The Mayer- Vietoris 
Lemma thus implies a G {h~^(A), h~^{A')). Since this is true for all r-perturbations h, 
we have a G U(r), as required. H 



Including intervals. We again need some definitions to unify the four cases into one. 
Given two intervals / and J of the same type, we say / includes into J, denoted as / ^ J, 
if f^^{I) includes as a pair in f^^{J). Unfolding the definition of the four types and 
assuming a < & < c < d, we have [5, c] ^ [a, d\, \b, d) ^ [a, c), (a, c] ^ (6, d], and 
(a, d) ^ {b, c); compare this with the Mayer- Vietoris diamonds in Figure|2] Suppose now 
that we have intervals I J, both of the same type. By the Interlevel Set Lemma, there 
are isomorphisms that take Hp(/~^(/)) and Hp(/~^( J)) to groups Gp{I) and Gp( J) in V. 
The inclusion induces a map on homology, which composes with these isomorphisms to 
give g : Gp(/) — > Gp( J). On the other hand, since the two groups are members of V, there 
is also a natural map from Gp(/) to Gp( J), namely the one that restricts to the identity on 
the span of their shared vectors and is zero otherwise. Not surprisingly, g is exactly that 
map. We formalize this claim and give a proof. 

Image Lemma. Let I ^ J and let Gp(/), Gp(J) be the corresponding p-dimensional 
groups in V. Then the image of g : Gp(/) — )• Gp( J) is a vector space in V, and its basis is 
in bijection with the multiset Wp{I) n Wp( J). 

Proof. To restate the lemma, we consider the diagram defined by the homology groups of 
the preimages of the including intervals, I ^ J, and the corresponding vector spaces in V: 

Hp(/-^(/)) ^ HpCrHJ)) 

t ; 

Gp(/) ^ Gp(J). 

The vertical maps are isomorphisms given by the Interlevel Set Lemma. The map h is 
induced by inclusion, and g maps a basis vector of Gp(/) to the same basis vector of 
Gp( J), if it exists, and to zero, otherwise. Hence, the basis of img consists of the vectors 
that are common to the bases of Gp(/) and Gp( J). This lemma states that we can get g by 
composing h with the two isomorphisms. Equivalently, the diagram commutes. To prove 
commutativity, we consider again the zigzag modules drawn as monotonic paths in the 
square; see Figured Since / J, we can find two non-crossing modules, one containing 
Hp(/~^(/)) and the other containing Hp(/^^(J)). To get a basis for imh, we translate 
intervals from one path to the other, keeping only the ones that cover both Hp(/~^(/)) and 
Hp(/~^ ( J)). Further translating these intervals to the hypotenuse gives the corresponding 
points in the persistence diagram. These points are precisely the ones shared by Wp(/) and 
Wp( J). In other words, im g in V is isomorphic to im h, as desired. H 
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Reading robustness. The Image Lemma allows us to compute the well groups and the 
well diagram associated to a single interval, / = (A, A'). The homology of f^^{I) can 
be read off the persistence diagram of /, as stated in the Interlevel Set Lemma. Simi- 
larly, the homology of (Xr,X^), where = ^rifh) and X^ = Xr(/^), can be read 
off the same diagram, as we now explain. By the Interval Formula, the well group for 
r is the intersection of the images of the maps t,, : Hp(Tr,T^) — > Hp(Xr,XJ,) and 
br : Hp(Br,B^) — > Hp(Xr,X^) induced by the inclusions. By the Image Lemma, this 
intersection corresponds to a pair of rectangles within the region of /^^(/); see the in- 
tersection between Wp (/) and the dotted rectangles in Figure |7] In the closed case, this 




H{f-\a,h]) m-Ha..b)) 



Figure 7: Reading the robust homology in the four different cases. The shaded region gives the basis 
of Hp(/^^(/)), while the dark shaded region gives the basis of the well subgroup, Up(r). 

intersection gradually recedes to infinity, while in the two half-open cases, the intersection 
disappears when r reaches half the length of the interval. Correspondingly, the well group 
shrinks gradually in the closed case, while it vanishes at or before r = (& — a)/2 in the 
half-open cases. Similarly, the well group vanishes when r reaches (6 — a)/2 in the open 
case. However, here it vanishes abruptly. More precisely, the range of the maps and b^, 
which is Hp(/~^(a + r,b — r)), approaches the homology group of the suspension of the 
level set at (a + b)/2, when r goes toward (6 — a) /2, before it suddenly becomes zero when 
r reaches that limit. 

In all four cases, a point contributes to the well group until r reaches a value at which 
the shrinking intersection no longer contains the point. Finding this value of r is easy since 
both rectangles shrink uniformly along all of their sides. Consider for example the case 
/ = [a, b] illustrated by the upper left design in Figure]?] For a point {x, y) G Dgm(/), the 
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value of r at which the point drops out of the relevant region is 



min{6 
min{6 
minja; 
min{2; 



x,y 
b, a 



a, a 



x,y 



b} if {x,y) eOi-d{f)nX[a,b], 

a} if {x,y) eExt{f)nX[a,b], 

y} if {x,y) eExt{f)n g[a,b], 

y} if {x, y) eRcl{f)n g[a,b]. 



The well diagram is the multiset of the values we get from the points in the persistence 
diagram. 

Measuring the difference. We can interpret the rank of the well group as a measure of 
the similarity between the image of the map : (T^, T^) — > (X^, X^) and the image of 
the map : (B^, BJ,) (X^, XJ,). Alternatively, we could use the cokemels of these two 
maps to measure their difference. Indeed, it is not difficult to prove counterparts of the 
Image Lemma for cokemels as well as for kernels. 

Co/kernel Lemma. Let I ^ J and let Gp(/), Gp(J) be the corresponding ^-dimen- 
sional groups in V. Then the kernel and cokernel of g : Gp(/) — !> Gp( J) are vector spaces 
in V, the basis of kcrg is in bijection with Wp(/) — Wp(J), and the basis of cokg is in 
bijection with Wp{J) - Wp{I). 

To measure the difference, we would therefore take the (algebraic) sum of the two coker- 
nels. Consider for example the open case. By the above lemma, we get a basis of coktr 
and cok br by setting J = (a + r, 6 — r) and first setting / to /i = {a + r,b + r) and second 
to /2 = {a — r,b — r). The basis of the sum, coktr + cok b^, is in bijection with the union 
of the two multisets of points, which is yVp( J) — Wp(/i) — Wpih)- 

5 Discussion 

The main contribution of this paper is the introduction of the point calculus for homol- 
ogy computations of level and interlevel sets. This comprises interlevel sets defined by 
closed, half-open, and open intervals, images, kernels, and cokernels of maps induced by 
inclusions, and the robustness of homology as defined by well groups. The point calculus 
provides a compact interface to a wealth of homological information that can be useful to 
researchers with and without background in algebraic topology. For the expert, it provides 
a compact summary of information that may be used to formulate conjectures about the 
topology of spaces and of functions. For the non-expert, the interface offers an intuitive ap- 
proach to understand the topology of datasets that by-passes the introduction of algebraic 
topology foundations. It is directly applicable to data in the form of continuous functions, 
which is common in medical imaging and in scientific visualization. 

We conclude by formulating an open question aimed at casting light on two- and higher- 
dimensional notions of robustness. This paper provides a solution to computing robustness 
when Y = M and perturbations are measured using the ioo-metric, and 13] shows that our 
results also hold for a broader class of metric function spaces. In |j9|, the authors give an 
algorithm when X is an orientable 2-manifold, Y = R^, and A is a point. Algorithms for 
other cases are not yet known. 
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